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By molecular dynamics simulation, the chemical vapor deposition of amorphous carbon onto graphite and 
diamond surfaces was studied. In particular, we investigated the effect of source H/C ratio, which is the ratio 
of the number of hydrogen atoms to the number of carbon atoms in a source gas, on the deposition process. 
In the present simulation, the following two source gas conditions were tested: one was that the source gas 
was injected as isolated carbon and hydrogen atoms, and the other was that the source gas was injected 
as hydrocarbon molecules. Under the former condition, we found that as the source H/C ratio increases, 
the deposition rate of carbon atoms decreases exponentially. This exponential decrease in the deposition 
rate with increasing source H/C ratio agrees with experimental data. However, under the latter molecular 
source condition, the deposition rate did not decrease exponentially because of a chemical reaction peculiar 
to the type of hydrocarbon in the source gas. 



1. Introduction 

Carbon is the most widely used material structure among all elements. It has two stable 
lattice structures, graphite and diamond. On the nanoscale, graphene, carbon nanotubes, 
carbon nanowalls, fullerene, and carbon onion are considered examples of nanostructures. 
The reason why carbon has many kinds of nanostructures is that sp 2 and sp 3 bonding states 
are balanced in binding energy. Silicon cannot form these kinds of nanostructure because its 
sp hybrid orbital consists of 3s and 3p orbitals, and then the sp 2 bonding state is weaker than 
the sp 3 bonding state. 

Amorphous structures are difficult to classify since exact quantities that can be used 
to characterize them have not yet been found. In particular, amorphous carbon is a complex 
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material because of the sp 2 and sp 3 bonding states. In addition, amorphous carbon mixed with 
hydrogen atoms is a more complex material. Thus, amorphous carbon has been investigated 
in numerous works. 1, 2 ) Chemical vapor deposition (CVD) and plasma CVD are often used in 
experiments. Diamond-like carbon (DLC) is considered as a useful material; thus, its formation 
process is actively investigated 3-6 ^ in engineering. However, the formation processes of the 
amorphous carbon and DLC by CVD are not known well. 

In the study of nuclear fusion, the amorphous carbon obtained by deposition, which is 
similar to the plasma CVD process, is disadvantageous. In nuclear fusion devices, hydro- 
gen plasma causes chemical erosion and chemical sputtering on divertor plates composed of 
carbon materials. Generated hydrocarbon molecules are transported by edge plasma, whose 
deposition creates amorphous carbon with hydrogen atoms. 7 ) The above phenomenon between 
plasma and surfaces is called plasma surface interaction (PSI). The PSI of the combination 
of plasma and carbon materials appears in various applications, such as the processing, coat- 
ing, and growth of thin films. We have investigated PSI on carbon materials with hydrogen 
atoms by molecular dynamics (MD) simulation with the clarification of nanoscale dynamics. 
In our previous works, the destruction mechanisms of the surface in PSI, especially chemical 
sputtering and chemical erosion, were investigated. 8 ' 9 ) In this study, we started to investigate 
the deposition process of a hydrocarbon gas onto carbon materials. On the nanoscale, the 
deposition in nuclear fusion devices and the CVD in laboratory experiments are considered 
to have similar mechanisms. 

Our simulation method and model are described in §2. In §3, we present and discuss the 
simulation results. This paper ends with conclusions in §4. 

2. Simulation Model 

2.1 Molecular dynamics for plasma surface interaction 

MD simulation for PSI is generally performed as follows. Figure 1 shows the conceptual 
diagram. A rectangular simulation box is considered. All particles move in the simulation 
box only. A target material is placed at the bottom of the simulation box. To make the 
surface of the target material vertical to the z-axis, the simulation box is placed relative to 
the coordination axes. The top and bottom of the simulation box are defined as positive 
and negative sides in the z coordinate. Incident particles are injected from the top side into 
the surface of the target material, where the z components of momentum of the injected 
particles are negative. The sizes of the simulation box in the x- and y-directions cannot be 
set arbitrarily because periodic boundary conditions in the x- and y-directions are imposed 
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on the simulation box. Therefore, the sizes of the simulation box in the x- and y-directions 
are determined as the lengths of the unit cells of integer pieces. If the surface slants toward 
the lattice direction, the unit cells are located in the simulation box at the slant. When the 
target material is selected, the species of atom composition, the lattice structure, and the 
Miller indices are determined. If the Miller indices are large, the size of the simulation box 
increases to line up the unit cells of integer pieces. The bottom of the target material should 
be fixed to the simulation box because the target material is entirely pushed by the incident 
particles. When the target material is a carbon material, flat (0001), armchair (1120), and 
zigzag (1010) surfaces of graphite, and diamond (100), (110), and (111) surfaces are typically 
used. The target material with an amorphous structure can be used after an amorphous 
material is created by the presimulation of the deposition or melting of a crystalline material. 

How physical phenomena can be treated strongly depends on the conditions of incident 
particles. If the incident energy is high, incident particles penetrate the target material deeply, 
which brings about the retention of particles, or they beat out particles in the surface area, 
which is called physical sputtering. Even if the incident energy is low, the covalent bonds in 
the surface area are broken by a chemical reaction with incident particles, which is called 
chemical sputtering or chemical erosion. It is difficult to classify the chemical sputtering and 
chemical erosion because the chemical sputtering does not beat out any atoms. If an incident 
particle includes an atom that can have many covalent bonds, e.g., a carbon atom, it creates 
deposits on the surface. 

2.2 Algorithm 

To deal with chemical reactions in MD, the modified Brenner reactive empirical bond order 
(REBO) potential was employed. 10 ~ 12 ) The modified Brenner REBO potential is given by 

U = E hS](^)-^(w,{^ B },{0)^](^)], 

i,j>i 

V minj) = + e W {-a [ij]nj ), (1) 

3 

v [tj) ( r ij ) = f[ij] ( r ij ) E B »M exp (~PnM Ti ) ' ( 2 ) 

n=l 

where is the distance between the i-th and j-tih atoms. The functions VSi an d Vj^i corre- 
spond to repulsion and attraction, respectively. The function by generates a multibody force. 
The modified Brenner REBO potential can treat structures peculiar to carbon, i.e., sp 2 and 
sp 3 bonds. The simulation time is developed using the second-order symplectic integration; 13 ) 
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the time step is 5.09 x 1CP 18 s. This time step is smaller than that of a general MD simula- 
tion. The modified Brenner REBO potential is represented by the interpolation between the 
potential of the bonding state and the potential of the unbonding state using cutoff functions, 
which have two cutoff lengths. Because the interpolation interval ranges from 0.3 to 0.6 A, a 
small time step is required. 

In this simulation, the material temperature was controlled using the Langevin equation: 

dU 

Pi = —Q^-lPi + Zi(t), (3) 

where and pi are the position and momentum of a carbon atom, respectively. The coefficient 
of friction 7 > 0. The random thermal force £j(i) satisfies 

m)) = 0, <&(<)&(f)> = 2D5 (t - *')<%, w 

where the symbol (• • •) indicates the expected value of a random variable, 5(t — t') is the Dirac 
delta function, and 5ij is the Kronecker delta. The magnification of the random thermal force 
D is given by the Einstein relation 

D = -ymik B T s , (5) 

where = 12 u is the mass of a carbon atom, k-Q is the Boltzmann constant, and T s is 
the setting temperature of a material. In the numerical simulation with the time step At, we 
substitute the difference equation 

\a r qjj / Ai\l 

^r(*)-w(*--2") At + ^ 2DAtB i (6) 

for eq. (3). The normally distributed random number Bi is generated by the Box-Muller 
transform 14 ^ from the uniformly distributed random number generated by single instruction 
multiple data (SIMD) oriented Fast Mersenne Twister. 15 ) In the present simulations, the 
thermal relaxation time I/7 was 0.1 ps. 

To investigate CVD by MD simulation, the simulation system was prepared as follows. 
The target materials were graphite and diamond. The surface of the graphite material was a 
flat (0001) surface of 2.00 x 2.17 nm 2 . The graphite material had a layer structure composed of 
four graphene sheets, which were stacked to form the "ABAB" pattern. Each graphene sheet 
consisted of 160 carbon atoms. The interlayer distance of graphite materials was set to 3.35 
A initially. The bottommost graphene sheet was regarded as the fixed layer, and then the 160 
carbon atoms in the fixed layer were bound by harmonic potentials whose centers were the 
initial positions of the atoms. Diamond materials with the (100) and (111) surfaces of 2.14 x 
2.14 and 2.02 x 2.19 nm 2 , respectively, were used. These diamond materials consisted of 648 
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and 640 carbon atoms, respectively. This difference in the number of carbon atoms made the 
sizes of the surfaces similar. In the fixed layers of these diamond materials, 72 and 80 carbon 
atoms were respectively bound by harmonic potentials. The graphite and diamond materials 
were connected to the Langevin thermostat. The kinetic energies of the atoms, which were the 
component atoms of the initial materials, were controlled by the thermostat, while those of 
the atoms injected into the surfaces were not controlled by the thermostat even if they were 
deposited on the surfaces. 

As the precursors of deposition, carbon atoms, hydrogen atoms, and hydrocarbon 
molecules were injected into the surfaces. The flux of the carbon atoms was 2.5 x 10 30 
atoms/m 2 s. The hydrogen atoms were always injected with the carbon atoms and the fluxes 
of the hydrogen atoms were from 2.5 x 10 29 to 2.0 x 10 31 atoms/m 2 s. The fluxes of hydrocar- 
bon molecules were set to have the same value in terms of the number of carbon atoms. For 
example, the fluxes of CH2 and C2H2 were 2.5 x 10 30 and 1.25 x 10 30 atoms/m 2 s, respectively. 
The incident energy was 1.0 eV. 

3. Results and Discussion 

First, only carbon atoms were injected into the diamond and graphite surfaces to compare 
the deposition mechanisms. All carbon deposits on the diamond and graphite surfaces grew 
with elapsed time linearly. The growth rates of the carbon deposits on the graphite (0001) 
surface and the diamond (100) and (111) surfaces were almost the same. The carbon deposits 
consisted of about 80% of sp 2 carbon atoms, 16% of sp 3 carbon atoms, and 5-8% of sp 1 
carbon atoms. Although the sp 3 carbon atoms were fewer than sp 2 carbon atoms, almost all 
the sp 3 carbon atoms were connected to each other (Fig. 2). The details of the compositions 
of the carbon deposits are shown in Table I. From the composition ratio, the carbon deposits 
produced by MD simulation were regarded as amorphous carbon (a-C). The carbon deposits 
were probably called graphite-like amorphous carbon, which is mainly composed of sp 2 carbon 
atoms. They were not diamond-like carbon (DLC), which has sp 3 carbon atoms of more than 
80%. 16 ~ 20 -' The density of the carbon deposits was about 2.4 g/cm 3 , which is higher than that 
of graphite but lower than that of DLC. These composition ratios and densities agreed with 
those of the amorphous carbon produced by tight binding MD. 21 ' 

To examine the dependence of the deposition mechanisms on the source H/C ratio, we 
performed MD simulation in which carbon and hydrogen were simultaneously injected as 
isolated atoms into the surfaces at a material temperature of 600 K. The ratio of incident 
flux of hydrogen atoms to that of carbon atoms was defined as the source H/C ratio. The 
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incident flux of carbon atoms was always fixed to 2.5 x 10 30 atoms/m 2 s, and the incident flux 
of hydrogen atoms was selected to be 0.1 to 8.0 times that of the carbon atoms. For example, 
when the flux of the hydrogen atoms was 1.25 x 10 30 atom/m 2 s, H/C = 0.5. In this paper, 
the deposition rate was defined as the ratio of the number of deposited carbon atoms to the 
number of injected carbon atoms. Figure 3 shows the deposition rate as a function of the source 
H/C ratio. The deposition rates plotted in this figure were calculated at the moment when 
2500 carbon atoms were injected. Because the surface area was about 4 nm 2 , the number of 
carbon atoms corresponding to a pure monolayer was 160. Then, a deposition rate less than 
0.064 (= 160/2500) indicates that the deposited carbon atoms could not form a sufficient 
amount of deposition layer, although the density of amorphous carbon is smaller than that 
of the pure crystal, strictly speaking. From this figure, we propose that the deposition rate of 
carbon atom Dq can be fitted by the exponential function of the source H/C ratio rjj/c as 



where the parameters (D , r ) were (0.698, 0.642), (0.663,0.6370), and (0.665, 0.441) for 
the diamond (100) and (111), and graphite (0001) surfaces, respectively. The parameter Ro 
was the deposition rate when only the carbon atom was the source atom. The deposition 
rates on the diamond (100) and (111) surfaces were almost the same. The deposition rate 
on the graphite (0001) surface was smaller than those of the diamond surfaces. In plasma 
CVD experiments, 22, 23 ** the relationship between deposition rate and source H/C ratio was 
researched. Actually, we confirmed that the deposition rate in the experimental data cited 
from Liu et al.[Fig. 3(c)], 22 ) which was defined as the deposition thickness per unit time, can 
be fitted by the exponential function of source H/C ratio. 

The compositions of the carbon deposits for the bonding state of carbon atoms were 
calculated as functions of the source H/C ratio, which are shown in Fig. 4. As the source H/C 
ratio increased, the number of sp 2 carbon atoms decreased but those of sp 1 and sp 3 carbon 
atoms increased slightly. In general, it is considered that sp 2 carbon atoms catch hydrogen 
atoms and change to sp 3 carbon atoms in a hydrogen-mixed layer. From the result, it is clear 
that hydrogen termination formed sp 1 carbon atoms. That is, sp 1 carbon atoms that have one 
C-H bond and one C-C bond were generated as carbon chain molecules, whose ends of carbon 
atoms are terminated by hydrogen atoms. The formation process from carbon and hydrogen 
atoms may prefer the use of the sp 1 carbon atom between C-C and C-H bonds. Actually, in 
hydrocarbon formation in the gas phase, carbon chain molecules are formed first. However, 



follows: 




(7) 
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Table I. Compositions of deposits. 



Surface 


sp 1 (%) 


sp 2 (%) 


sp 3 (%) 


Density (g/cm 3 ) 


Diamond (100) 


6.5 


77.3 


16.2 


2.49 


Diamond (111) 


7.5 


75.8 


16.7 


2.34 


Graphite (0001) 


4.8 


78.7 


16.5 


2.42 



after carbon chain molecules are formed in the gas phase, they change into ring and cluster 
structures slowly. Prom this fact about the formation process in the gas phase, if the present 
deposition simulation is executed for a long time at a lower incident flux, sp 1 carbon atoms 
can potentially change into sp 2 or sp 3 carbon atoms. 

In the deposition of the amorphous carbon material, the source gas not fully resolved 
into atoms, the deposited particle could have possibly been molecules or radicals. Then, the 
deposition using C2, CH, CH2, C2H2, and carbon chain C4H2 molecules was tested. Figure 5 
shows the deposition rates of carbon atoms in the cases of these hydrocarbon molecules. The 
deposition rate of C2 was greater than that of the carbon atom, while the deposition rates 
of the hydrocarbon molecules were smaller. C2H2 hardly formed carbon deposits on all the 
surfaces. For CH and CH2, the deposition rates on the diamond (100) surface were greater 
than that on the diamond (111) surface. Only on the graphite (0001) surface, the deposition 
amount of CH2 was passably larger than that of CH. Here, the reason why CH2 could deposit 
more easily than CH although the source H/C ratio of CH2 was greater than that of CH was 
determined by comparing CH2 deposition with CH deposition. Figure 6 shows the deposition 
process actually observed in the present MD simulation. The photograph of CH injection in 
Fig. 6(a) illustrates many C2H residues, whose two carbon atoms are in the spi state, on the 
graphite surface. A typical CH deposition process is as follows [Fig. 6(c)]. First, a CH molecule 
was adsorbed on the surface and formed a CH residue. Next, the second CH molecule reacted 
with the CH residue and then a C2H residue was generated. The C2H residue could not grow 
because the hydrogen atom on top of the C2H residue blocked the third CH molecule. On 
the other hand, the photograph of CH2 injection in Fig. 6(b) shows many sp 3 carbon atoms 
that have two C-C bonds and two C-H bonds. In the CH2 deposition process [Fig. 6(d)], CH2 
molecules form a polymer-like structure on the graphite surface. Because the polymer-like 
structure had alternating trans and cis isomers, the structure was not straight. That is, the 
hydrogen atom, which reflected incoming source particles, was not located on only the top of 
the polymer. Therefore, CH2 could be adsorbed on top of the polymer one after another. 
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4. Conclusions 

The chemical vapor deposition of amorphous carbon was investigated by MD simulation. 
The base surfaces of the present simulation were the graphite (0001) and diamond (100) and 
(111) surfaces, and the source gas was injected in the following two states: isolated carbon 
and hydrogen atoms being injected according to the source H/C ratios of 0.0 to 8.0, and 
hydrocarbon molecules being injected as CH, CH2, C2H2, C4H2, and C2. We determined the 
dependence of the source H/C ratio on the deposition process of amorphous carbon. As a 
result, in the isolated atom injection, we found that the deposition rate of carbon atoms 
can be fitted by the exponential function of the source H/C ratio. It is confirmed that the 
experimental data cited from Liu et al. 22 -* also agree with this exponential decrease in the 
deposition rate with increasing source H/C ratio. However, in the hydrocarbon molecular 
injection, the deposition rate exponentially decreases as the source H/C ratio increases. As 
an example of the breaking of the exponential decrease, we chose the MD simulation result 
showing that CH2 was more rapidly deposited than CH although the H/C ratio of CH2 was 
greater than that of CH. We expect that the exponential decrease in the deposition rate 
with increasing source H/C ratio tends toward macroscopic modeling for CVD of amorphous 
carbon. 
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Figure Captions 

Fig. 1. (Color online) MD simulation system for PSI. 

Fig. 2. (Color online) Photograph of the MD simulation on the diamond (100) surface. The 
yellow, green, and red spheres indicate the sp 1 , ap 2 , and sp 3 carbon atoms, respectively. 
The blue spheres indicate carbon atoms that have one covalent bond. The purple spheres 
indicate carbon atoms connected by five covalent bonds or more. 

Fig. 3. (a) Deposition rates of the carbon atoms on the diamond (100) and (111), and graphite 
(0001) surfaces as functions of the source H/C ratio and (b) its enlargement. Data points 
under the dashed line of 0.064 in (b) are depositions less than the number of carbon atoms 
corresponding to a monolayer, (c) Experimental data cited from Liu et al. 22 ) 

Fig. 4. The compositions of the carbon deposits for the bonding state of carbon atoms were 
calculated as functions of the source H/C ratio: (a) a diamond (100) surface, (b) a diamond 
(111) surface, and (c) a graphite (0001) surface. 

Fig. 5. Deposition rates of carbon atoms in the cases of the hydrocarbon molecules: (a) a 
diamond (100) surface, (b) a diamond (111) surface, and (c) a graphite (0001) surface. 

Fig. 6. (Color online) Photographs of MD simulation in the cases of (a) CH and (b) CH2 
injections. The blue, yellow, green, and red spheres indicate atoms having one to four 
covalent bonds, respectively. Typical deposition processes of (c) CH and (d) CH2. The 
white and gray spheres indicate hydrogen and carbon atoms, respectively. 
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